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^ ; Abstract 

In a recent paper Dar and Shaviv (1994) presented results of calculations 
with an improved standard solar model (SSM) that suggest a standard 
physics solution to the solar neutrino problem. In a subsequent publi- 
cation with a similar title Bahcall et al. (1994) have claimed to refute 
the results obtained by Dar and Shaviv. Here we show that the crit- 
icism and conclusions of Bahcall et al. are unjustified. In particular, 
0\ '. their attempts to reproduce the results of Dar and Shaviv failed because 

■ they did not use the same nuclear cross sections and did not include in 
their calculations major physics improvements which were included in 
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the improved SSM code used by Dar and Shaviv. 



Introduction 

To facilitate understanding our response to a critique on our paper "A 
Standard Physics Solution To The Solar Neutrino Problem ?" (Dar and 
Shaviv, Phys. Rev. Lett. Submitted 1994), posted by 15 authors! 
(Bahcall et al. 1994) on this electronic bulletin board, we quote here (in 
italics) the main comments and criticism of these authors and insert our 
responses (in Roman fonts) between them. 

Bahcall et al. abstract : 



"The claim by Dar and Shaviv that they have found a standard model 
solution to the solar neutrino problem is based upon an incorrect assump- 
tion made in extrapolating nuclear cross sections and the selective use 
of a small fraction of the nuclear physics and of the neutrino data. In 
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addition, five different solar model codes show that the rate obtained for 
the chlorine experiment using the Dar-Shaviv stated parameters differs 
by at least 14<j from the observed rate. " 

DS Response: 

Dar and Shaviv did not claim to have solved the CI solar neutrino prob- 
lem. Rather, they claimed agreement, within theoretical and ex- 
perimental uncertainties, between the predictions of an improved 
standard solar model and all solar neutrino observations after 1986 
(i.e., with all the published results of Kamiokande, SAGE and GALLEX 
and with the Homestake observations after 1986). This fact strongly 
suggests a standard physics solution to the solar neutrino problem. 

Dar and Shaviv did consider all relevant nuclear data and did not intro- 
duce incorrect assumption in extrapolating nuclear cross sections to low 
energies (see below). 

The five different SSM codes quoted by Bahcall et al. do not and should 
not reproduce the results of Dar and Shaviv when using their derived 
reaction rates for only two nuclear reactions and their value for the so- 
lar luminosity (Particle Data Group 1992): Although Dar and Shaviv 
have stated explicitly in their paper that for all the reactions, except for 
7 Be(p,7) 8 B and 3 He(a,7) 7 Be, they have used the reactions rates com- 
piled by Caughlan and Fowler (1988), Bahcall et al. chose to "reproduce" 
the results of Dar and Shaviv with different reaction rates. Moreover, 
the Bahcall et al. codes include incorrect physical effects and unjustified 
physical approximations which were not used in the Dar-Shaviv calcula- 
tions. In particular, these codes include incorrect screening enhancement 
of nuclear reaction rates in the Sun (omitted in the Dar-Shaviv code), 
assume nuclear equilibrium for the CNO bi cycle (not assumed in the 
Dar-Shaviv code), impose complete ionization of all elements everywhere 
in the Sun in calculating diffusion and the equation of state (not imposed 
in the Dar-Shaviv SSM code), replace a complex variety of heavy ele- 
ments by a single effective element in calculating diffusion, opacities and 
the equation of state (80 different isotopes are included explicitly in the 
Dar-Shaviv SSM code. In particular, the Dar-Shaviv code follows the 
destruction of the light elements from Pre Main Sequence composition 
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till the present day) and do not calculate the evolution of the Sun during 
the Pre Main Sequence stage (the Dar-Shaviv SSM calculations include 
the evolution of the Sun during the Hayashi phase). Most surpris- 
ingly, Bahcall et al. have ignored the fact that the Bahcall- 
Pinsonneault results were reproduced quite accurately by the 
code used by Dar and Shaviv when similar physical assump- 
tions and approximations were used, i.e. before introducing 
the above improvements (see Kovetz and Shaviv 1994). 

In their paper Dar and Shaviv have stated explicitly that the results 
of the Homestake experiment prior to 1986 are inconsistent with their 
improved SSM predictions but they also pointed out that these results 
look statistically incompatible with the Homestake results after 1986, 
assuming a constant Sun. The 14a discrepancy, between the Bahcall et 
al. "reproduction" of the Dar-Shaviv prediction and the observed rate 
in the CI experiment averaged over 23 years, is partly due to improper 
reproduction of the Dar-Shaviv results by Bahcall et al. and partly due 
to not including in a the large theoretical uncertainty and the most 
probably large unknown systematic errors in the CI results (as implied 
by the strong variation of the results as function of time). 

The Homestake results after 1986 are compatable with those obtained by 
Kamiokande during the same time period. Furthermore, the results of 
GALLEX and SAGE can be predicted correctly from these observations 
and the solar luminosity using essentially only standard physics conser- 
vation laws (shown in the Dar-Shaviv paper). Therefore, Dar and Shaviv 
concluded that the solar neutrino experiments do not provide evidence 
for physics beyond the standard electroweak model and that only future 
experiments like SNO and Superkamiokande will be able to rule out a 
standard physics solution to the solar neutrino problem. 

Bahcall et al. claim : 

"2. Choosing Part of the Chlorine Solar- Neutrino Experimental Data 

Dar and Shaviv chose to consider (see their Figure 1) only four years 
of data from the chlorine solar neutrino experiment beginning in 1987, 
although 23 years of data have been reported. For the period beginning in 
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1987, the experimental measurement of the chlorine rate is 2.8±0.3 SNU. 
The measured rate for the entire period for which data has been reported 
is 2.28±0.23 SNU. The measured rate in the chlorine experiment, during 
the short period considered by Dar and Shaviv, is still more than 4<j less 
than their calculated result of 4.2 SNU and is more than 14<r from the 
value we calculate for their stated parameters. " 

DS Response: 

Dar and Shaviv stated clearly in their paper (in the abstract, text and 
conclusions) that they compared their predictions with all the reported 
solar neutrino observations after 1986, i.e., with all the reported results 
from Kamiokande, SAGE and GALLEX and from runs 90-120 of the 
Homestake CI experiment during the time period 1986.8-1992.4, and not 
with "only four years of data" . In Fig. 1 they showed the data published 
by Homestake and Kamiokande for only the period 1987-1990 in order 
to avoid showing experimental data that was privately communicated to 
them by the authors prior to their publication. 

The 37 Ar average production rate observed in the CI experiment in 
runs 90-120 (during the time period 1986.8-1992.4), which has been an- 
nounced in recent international conferences and has been communicated 
to us also by K. Lande (who is one of the co-authors of the paper by 
Bahcall et al.), is 0.61 ± 0.5 atom/day. Since 1 atom/day = 5.35 SNU, 
this rate corresponds to 3.26 ±0.28 SNU before background subtraction 
and to 3.01 ± 0.28 SNU after background subtraction of 0.047 ± 0.016 
atom/day. The last number is the present best estimate of the back- 
ground in the CI experiment, reported by Ray Davis at two recent in- 
ternational meetings: "Neutrino Telescopes 1994" in Venice and "The 
Solar Neutrino Problem - Astrophysical or Particle Physics Solution" in 
Grand Sasso. 

The average production rate of 37 Ar in CI during the time period 1970.8 
- 1985.3 observed by the Homestake experiment, as was communicated 
to us by K. Lande, is 0.47±0.04 atoms/day which corresponds to 2.09 ± 
0.27SNU after 0.08±0.03 atoms/day background subtraction, or to 2.26± 
0.023 SNU after 0.047 ± 0.016 atoms/day background subtraction. 

As explained in our paper, the Homestake results prior and after 1986 
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look statistically incompatible with a steady state Sun. The annual 37 Ar 
production rate has increased by more than 50% since the beginning of 
the Homestake experiment. One can ignore that and compare theory 
with the rate averaged over the whole live time of the Homestake exper- 
iment, assuming that the difference between the results before and after 
1986 is due to a statistical fluctuation (with a rather small probability) 
as advocated by Bahcall et al. and by some other authors. Our ap- 
proach is different. Kamiokande began its Solar Neutrino observations 
in 1987. It is only natural to compare the results of the Homestake and 
Kamiokande experiments during the time period when both were run- 
ning, in particular, in view of the fact that Ray Davis and other authors 
proposed that solar neutrino flux is time dependent and anticorrelates 
with the solar activity. It is worth to note, and perhaps even significant, 
that after the installment of the new pumps in the CI experiment in 1986 
the results from the new CI experiment (Homestake II, runs 90-120 dur- 
ing 1986.8-1992.4) are not in conflict with the results from Kamiokande. 
Moreover, during this time period neither the Homestake results nor the 
Kamiokande results showed any statistically significant anticorrelation 
with solar activity, nor any significant time variation, which have been 
claimed with high significance for the CI data prior to 1986 by a few 
authors of the Bahcall et al. paper and by other authors. 

In view of all these facts and in view of the agreement between our SSM 
predictions and the Kamiokande results for the 8 B solar neutrinos, it is 
scientifically unjustified to overlook the possibility that the results after 
1986 are closer to the truth than both (a) the results before 1986 and 
(b) the weighted average of the results before and after 1986. 

The Dar-Shaviv prediction for the CI experiment is 4.2 ± 0.6 SNU. The 
results from the uncalibrated CI experiment after 1986 correspond to 
3.01 ±0.28 SNU, where the error is statistical only, not including possible 
systematic errors. Bahcall et al. concluded that the experimental result 
is 4<7 below the theoretical prediction ignoring the theoretical uncertainty 
and possible systematic errors. We, however, consider the above results 
to be consistent within the theoretical and experimental uncertainties. 

The impressive 14<r discrepancy between the expected rate in the CI ex- 
periment, as calculated by Bahcall et al. with the Dar-Shaviv nuclear 
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parameters, and the average rate at Homestake obtained for the whole 
23 years of observations is partly due to the incorrect physical approxi- 
mations used in the Bahcall et al. calculations (see also below), partly 
due to ignoring the large uncertainty in the theoretical prediction, and 
perhaps due to possible unknown large systematic errors in the uncali- 
brated CI experiment (as implied by the significant increase of the rate 
as a function of time since the starting of the CI experiment). 

Bahcall et al. claim. 

"3. Extrapolating Nuclear Cross Sections with an ad hoc Assumption 
and Using Only Part of the Experimental Data 

The theoretical models that have been used previously Johnson et al. 
1992, Kajino et al. 1984, Parker and Rolfs 1991 to extrapolate the 
cross section data to solar energies properly take account of the finite 

nuclear size effects along with the other effects These effects include 

nuclear structure, the strong interaction, energy dependent operators in 
the transition matrix elements, antisymmetrization between the colliding 
nucleons, finite nuclear size, the final-state phase space, and the contri- 
butions from other partial waves. . . . 

... Dar and Shaviv assumed incorrectly that the only energy- dependent 
effect besides point-nuclei barrier penetration is nuclear size... Dar and 
Shaviv did not take account of all the energy dependencies nor of all 
the available nuclear physics data... Their selective use of data and this 
incorrect assumption explain why the Dar and Shaviv answers for S(0) 
differs from the standard values obtained by nuclear physicists ..." 

DS Response: 

Bahcall et al. describe a wishful idealistic situation, where nuclear reac- 
tion theory is an exact and tested theory which correctly predicts all cross 
sections (or at least their energy dependence over a wide range). This sit- 
uation is far from reality. Actually, nuclear reaction theory is an approx- 
imate theory which generally provides only a reasonable parametrization 
of measured nuclear cross sections over a limited range of energy using 
many free parameters that are directly adjusted to fit the data. In fact, 
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for most reactions, including the major reactions in the Sun, it cannot 
even distinguish which of the different measurements of a cross section 
is the correct one (normalization as well as energy dependence). To 
make our point more concrete, let us consider the model of Johnson et 
al. (1992)) used "by the nuclear physicists" (together with an ad hoc 
procedure) to extrapolate the various measurements of the cross section 
for 7 Be(p,7) 8 B to zero energy to obtain S17(0)=22.4 eV-b, the value 
advocated by Bahcall et al. In spite of all the good features of that 
model, which were listed by Bahcall et al., Johnson et al. (1992) have 
not demonstrated that their model correctly predicts: 

(i) the p-wave resonance (i.e. its position, width and magnitude), 

(ii) the magnitude of the cross section, 

(iii) the energy dependence of the cross section over the entire energy 
range where it has been measured. 

In fact Johnson et al. did not trust their model energy dependence 
(which they refer to as "suspect theoretical calculations above 430 KeV" 
cf. the last paragraph in Ap.J. 392,320,(1992) page 325) to extrapolate 
the cross sections measured above the resonance energy by Kavanagh 
1960, Parker 1966, 1968 and Vaughn et al. 1970, and choose instead 
to extrapolate these data to low energies below the resonance according 
to the energy dependence of the averaged Kavanagh et al. 1969 and 
Filiponne et al. 1983 experiments (an ad hoc prescription). 

Moreover, 

(iv) for the 7 Li(n,7) 8 Li capture cross section Johnson et al. pre- 
dicted 30.6 mb while the experimental values are 40.2 ± 2 mb (Imhof et 
al. 1959) and 45 ± 3.0 mb (Lynn et al. 1991), 

(v) for the 1=2 7 Li+n isotopic analog channel Johnson et al. pre- 
dicted a scattering length a 2 = 0.26 /m, while the experimental value is 
a 2 = -3.59 ±0.06 fm . 

This situation is rather typical (see for instance Descouvemont and 
Baye 1994, in particular Figs. 5.6. therein). 

The above points are not intended to criticize the otherwise good works 
of Johnson et al., and Descouvemont and Bay, but rather to emphasize 
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our contention that since no exact theory exists for either the compli- 
cated quantum mechanical nuclear many body problem or the direct 
nuclear reactions, one should try to rely as much as possible on "model 
independent" features rather than on the choice of a specific model for 
extrapolating cross sections to very low energies. 

Three such features were used in our paper: 

(1) The measured position, width and height of a low energy reso- 
nance in a given partial wave together with the effective range approxi- 
mation determine uniquely the contribution of that partial wave to the 
total cross section all the way down to zero energy (Kim et al., 1994). 
It was used by us to subtract the p-wave contribution from the cross 
section for 7 Be(p,7) 8 B. 

(2) Because of the Coulomb and centrifugal barriers, the relative 
contribution of different non resonating partial waves to the total cross 
section, at energies well below the Coulomb barrier, are almost model 
independent, while the magnitude and energy dependence of the indi- 
vidual partial waves are model dependent. Thus, we used only the ratios 
calculated by various groups (but neither the magnitude nor the energy 
dependence) to extract from the measured cross section the s-wave con- 
tribution. 

(3) Bound state wave functions decrease exponentially outside the 
nucleus while the incident Coulomb wave functions for energies well be- 
low the Coulomb barrier decrease exponentially from the classical turn- 
ing point towards the nuclear surface. Consequently, most of the con- 
tribution to the overlap integrals comes from the vicinity of the nuclear 
surface. The s-wave cross section therefore, is proportional to the abso- 
lute magnitude squared of the incident Coulomb wave function near the 
nuclear surface, rather than at the origin, i.e. to the Gamow (WKB) 
barrier penetration factor for finite nuclear radii rather than the Gamow 
barrier penetration factor for point nuclei. 

When extracting this energy dependence from the s-wave contributions 
to the measured cross sections for the three reactions, 7 Be(p,7) 8 B, 3 He(ai, 7) 7 Be| 
and 3 He( 3 He,2p) 4 He, we found no evidence for any significant energy de- 
pendence of our defined S(E) for E well below the Coulomb barrier. The 
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absence of any noticeable energy dependence in S(E) extracted from the 
measured cross sections for the above three reactions, while S(E) shows 
strong energy dependence, which is not completely accounted for by the 
complicated theoretical models, is in our opinion, a further support for 
our simple extrapolation procedure. 

Bahcall et al. claim : 

"For the determination of the low- energy cross section factor, S34, for 
the 3 He(a, / y) 7 Be reaction, Dar and Shaviv apparently adjusted the ra- 
dius parameter R so that the energy dependence of S(E) is mostly re- 
moved for two of the nine (Parker and Rolfs 1991) existing experiments. 
(They seem not to have noticed that the value of R that they obtain is 
very different from the measured radius of 2.8 fm determined by elec- 
tron scattering.) They did not allow for the other energy dependencies 
discussed above and they only took account of two of the experiments. " 

DS Response : 

Dar and Shaviv did not adjust the radius parameter R = i?( 3 He)+i?( 4 He)| 

so "that the energy dependence of S(E) is mostly removed for two of the 
nine existing experiments". Dar and Shaviv used R = _R( 3 He)+i?( 4 He)~ 
3.6 fm where i?( 3 He)~ 2.0 fm and _R( 4 He)~ 1.6 fm were determined 
from electron scattering measurements (e.g., Strueve et al. 1992; Wu et 
al.1994) and from scattering of high energy strongly interacting particles 
(e.g. Tanihata et al. 1988), respectively, to calculate S3j(E) for all nine 
experiments (the world data). The straight line in Fig. 2 represents the 
weighted average of S37 (E) for these world data. (Dar and Shaviv were 
aware that R( 7 Be) ~ 2.8 fm, but that is not the appropriate radius to 
be used). 

Bahcall et al. claim : 

"Dar and Shaviv cited the preliminary Coulomb dissociation work de- 
scribed in preprint form (Motobayashi et al. 1994) as evidence for a 
lower- than- standard value for the crucial cross section factor for the 
7 Be(p, / ~f) 8 B reaction. When the E2 contribution to this reaction is 
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taken into account (Langanke and Shoppa 94), the preliminary Coulomb- 
dissociation value differs from the six direct measurements (Parker and 
Rolfs 1991) of the 7 Be(p,"f) 8 B cross section by a factor of two while the 
estimated uncertainty (Johnson et al. 1992) in direct measurements is 
only 11%. Moreover, there are still some unanswered questions about 
the application of the Coulomb -dissociation method for determining ra- 
diative capture cross sections, aside from the experimental difficulties 
inherent in covering a sufficient range in energy and angle to validate 
the reliability of any inferences. " 

DS Response: 

The results described in the paper by Motobayashi et al. (1994) that 
was submitted for publication in Physical Review Letters and was cited 
by Dar and Shaviv were shown (Fig. 3 of Dar and Shaviv ) to be in 
good agreement with the "world average" s-wave cross section that was 
extracted from the four measurements of the 7 Be(p, 7) 8 -B reaction in the 
same energy range, (which is smaller there by approximately a factor two 
than the total cross section) and in particular with that extracted from 
Vaughn (1970) and Filippone (1983) if one allows a small (< 10%) Ml 
contribution. The world data that includes that of Motobayashi et al. ( 
1994) yield 517(0) = 17 ±2 eV b. The data of Motobayashi et al. yields 
approximately the same value. 

Neither the suggestion that the results of Motobayashi et al. (1994) are 
preliminary nor their analysis by Langanke and Shoppa (1994) justify 
their omission from the "world data" that is used to determine SI 7(E). In 
fact, the analysis of Langanke and Shoppa of the Motobayashi et al. data 
is highly questionable. They used as input E2 nuclear matrix elements 
calculated by others (Kim, Park and Kim, 1987). But if instead one uses 
the E2 nuclear matrix element calculated by Descouvemenot and Baye 
(1994) one gets an effect that is approximately 2.4 times smaller at 0.6 
MeV. In addition, Langanke and Shoppa find the largest E2 contribution 
at 0.6 MeV, which is on the Ml resonance. The E2 contribution of 
the resonance is of no direct relevance to the low energy (s-wave) cross 
section. 

In the data points of Motobayashi et al. (1994), measured at 0.8 
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and 1.0 MeV, that Langanke and Shoppa analyzed, there is essentially 
no evidence for an E2 component, as they state themselves. In fact, in 
all 15 or so data points shown by Motobayashi et al. (1994) there is only 
one point at approximately 4°, 0.6 MeV bin (on resonance) where the 
data deviates by approximately 2 a from a predicted pure El behavior 
(the E2 angular dependence is sufficiently different from the El's). A 
claim based on that one point (out of 15 shown and some 30 measured in 
total [Motobayashi et at. private communication] that there is evidence 
for E2 contribution in the Motobayashi et al. data is a bit too much. 

Bahcall et al. claim : 

"For other nuclear reactions, Bar and Shaviv have used the low- energy 
cross section factors from an earlier review ( Caughlan and Fowler 1988) 
which provided fitting formulae suitable for use at temperatures (~ 10 9 
K) much higher than are reached in the Sun (~ 10 7 K). The quantita- 
tive effect of these approximations is difficult to estimate, especially since 
other authors ( see section 4 ) use-for solar calculations- explicit formulae 
that are suitable for the lower solar temperatures. However, an approxi- 
mate discussion of using the fitting formulae for higher temperatures at 
solar temperatures has been given (Bahcall 1992); the principal effects of 
the high-temperature formulae are in the direction to decrease the pre- 
dicted 7 Be and 8 B neutrino fluxes." 

DS Response: 

We have verified that for solar temperatures the reaction rates for all the 
reactions, except 7 Be(p,7) 8 B and 3 He(oi, 7) 7 Be, which were compiled by 
Caughlan and Fowler (1988) are not significantly different from those 
obtained by Bahcall and Ulrich (1989) and by us. In fact the Kovetz- 
Shaviv SSM code with the Caughlan - Fowler (1988) nuclear reaction 
rates (for all reactions) and all other assumptions and input as used by 
Bahcall and Pinsonneault 1992 reproduces their predicted pp, pep, hep, 
7 Be and 8 B solar neutrino fluxes within 2% accuracy (see column KS of 
Table I of the Dar-Shaviv paper or Kovetz and Shaviv 1994). 

Bahcall et al. comment : 
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"Dar and Shaviv discuss at some length the fact that the Debye-Huckel 
approximation to the screened nuclear potential is not correct everywhere 
in the Sun. It is not clear what they recommend (although they say the 
effect is small), nor if they are aware that modern screening calculations 
go well beyond what they discuss..." 

DS Response: 

Dar and Shaviv discuss at some length the reasons why the screening 
enhancement of the nuclear reaction rates near the center of the Sun 
are negligible. The conventional screening enhancement factors used by 
Bahcall et al. in their SSM codes should be taken out. This modifi- 
cation in the Bahcall-Pinsonneault SSM code will reduce significantly 
their predicted 8 B solar neutrino flux (about 15%). It does not change 
significantly the predictions if this modification is introduced after im- 
plementing all our other suggested improvements. 

Bahcall et al. 

4- Solar Model Calculations 

Six authors of this paper (Bahcall, Christens en- Dalsgaard, DeglTnnocenti,^ 

Glasner, Pinsonneault, and Proffitt) have repeated the solar model cal- 
culations of Dar and Shaviv using the non-standard parameters that Dar 
and Shaviv chose, namely, a solar luminosity of 3.826 x 10 33 erg s _1 
and low energy cross-section factors of £34(0) = 0.45 KeV-b and 
#17(0) = 17 eV-b. Dar and Shaviv did not specify in their preprint 
many of the important input quantities in their model; they did not state 
what they used for the element abundances, the radiative opacities, the 
equation of state, and the neutrino cross sections. They did not say 
which of the several available prescriptions for diffusion they used. We 
have therefore carried out calculations using a variety of different choices 
for these quantities, namely, the choices made previously as their best 
estimates by the six different authors who used five independent stellar 
evolution codes ... None of the well-tested solar codes that we have used 
are able to reproduce the Dar and Shaviv results for neutrino fluxes. 

DS Response: 
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The input physics used by Dar and Shaviv is the most updated one. The 
detailed SSM code with all input physics which was used by Dar and Sha- 
viv is documented in detail in the paper of Kovetz and Shaviv (Ap. J, 
May 1, 1994) which is cited by Dar and Shaviv. The Kovetz-Shaviv 
paper is and was available in a preprint form upon request. Bahcall et 
al. apparently overlooked this paper and many physical improvements 
introduced by Dar and Shaviv in the SSM calculations. All the five dif- 
ferent SSM codes of Bahcall et al. do not and should not reproduce the 
results of Dar and Shaviv when using the same measured solar lumi- 
nosity that was used by Dar and Shaviv and the same nuclear reaction 
rates for two reactions only, 7 Be(p,7) 8 B and 3 He(a,7) 7 Be, because 
these SSM codes use reaction rates different from the Caughlan-Fowler 
(1988) reaction rates used by Dar and Shaviv for all other reactions, 
and because they include incorrect physical effects and unjustified ap- 
proximations which were avoided in the Dar-Shaviv calculations: They 
include incorrect screening enhancements of nuclear reaction rates in 
the Sun (which were not included in the Dar-Shaviv calculations), as- 
sume (some of them) nuclear equilibrium for the CNO bi cycle (which 
was not assumed by Dar and Shaviv), impose complete ionization of all 
elements everywhere in the Sun in calculating diffusion and the equa- 
tion of state (which was not imposed by Dar-Shaviv), replace a complex 
variety of heavy elements by a single effective element in calculating 
diffusion, opacities and equations of state (while Dar-Shaviv considered 
separately each of the 80 isotopes) and do not calculate the evolution of 
the Sun during the pre main sequence stage (while Dar and Shaviv did). 
Most surprisingly, Bahcall et al. have ignored the fact that the Bahcall- 
Pinsonneault results are reproduced quite accurately by the code used by 
Dar and Shaviv when they imposed the same physical approximations 
that were used by Bahcall and Pinsonneault. (The Dar-Shaviv code 
also includes several improvements and enhancements in the accuracy of 
evolutionary models calculations never applied before). 

In particular, Bahcall (1992) found the following empirical depen- 
dence of the 8 B solar neutrino flux on the cross sections for the major 
nuclear reactions in the Sun: 

A /8t>\ -2.6 -0.4 n.8 1.0 

0f©( ~ ^11 ^33 a 34 <?17 ■ 
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Thus, the use of our values (MeV-b units), Su = 4.07 x 10 22 , £33 = 
5.60 x 10 3 S 34 = 0.45 and S 17 = 0.017 instead of the values Su = 
4.00 x 10" 22 , S33 = 5.0 x 10 3 , S34 = 0.533 and S 17 = 0.0224 which were 
used by Bahcall et al. should decrease the Bahcall - Pinsonneault (1992) 
prediction for 4> UQ ( 8 B) by a factor « 0.61 . Similarly, the omission of 
the weak screening enhancement factors, exp(U(0)/kT) where U(0) ~ 
ZiZ 2 e 2 /Rd-, from the nuclear reaction rates should decrease the 8 B neu- 
trino production rate near the center of the Sun (T c « 1.571 x 10 7 K 
and Rd ~ 2.8 x 10 _9 cm) by a factor « 0.89 . Thus, the combined effect 
of the use of the Dar-Shaviv reaction rates and the omission of the weak 
screening enhancement, according to Bahcall 1992, is to reduce the Bah- 
call - Pinsonneault (1992) prediction, cf) UQ ( 8 B) « 5.69 x 10 6 cm" 2 s _1 , 
by a factor w 0.54, i.e., to </v ( 8 B) « 3.09 x 10 6 cm s . The reduced 
solar luminosity brings it down by additional 5% and all other improve- 
ments/changes introduced by Dar and Shaviv should bring it further 
down (to their calculated value, l/0 ( 8 B) « 2.77 x 10 6 cm" 2 s _1 ?). 

In view of all these, we find no point in discussing the conclusions of 
Bahcall et al. which are based on misleading comparisons. Bahcall et 
al. should have introduced all the suggested improvements/changes in 
their independent codes and then try to reproduce the Dar-Shaviv results 
before drawing their conclusions. 

Bahcall et al. 

5. Conclusions: 

" Dar and Shaviv did not succeed in solving the solar neutrino problem... 
Their failure to solve the problem is not surprising since it has been 
demonstrated elsewhere (Bahcall and Bethe 1993) that the chlorine and 
the Kamiokande experiments are inconsistent another (if one assumes 
standard electroweak theory)..." 

DS Response: 

Dar and Shaviv did not claim to have solved the CI solar neutrino 
problem. They claimed agreement, within theoretical and experimen- 
tal uncertainties, between the predictions of an improved standard solar 
model and all solar neutrino observations after 1986 (i.e., with all the 
published results of Kamiokande, SAGE and GALLEX and with the 
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Homestake observations after 1986). That suggest a standard physics 
solution to the solar neutrino problem. Moreover, the Homestake results 
after 1986 are consistent with those obtained by Kamiokande during the 
same time period. Their results together with the observed solar lumi- 
nosity predict correctly the observed production rate in GALLEX and 
SAGE, using essentially only standard physics conservation laws (see 
Dar and Shaviv 1994). Finally, the spectrum of 8 B neutrinos measured 
by Kamiokande is not different from that predicted by the standard elec- 
troweak model. In view of all these facts Dar and Shaviv concluded that 
the solar neutrino experiments do not provide solid evidence for physics 
beyond the standard electroweak model. Only future experiments with 
large statistics and precise measurements of the energy spectrum and the 
lepton flavor content of solar neutrinos, such as the Super kamiokande 
light water experiment and the SNO heavy water experiment will be 
able to provide reliable evidence for physics beyond the standard elec- 
troweak model, which has not been found so far in all precision tests of 
the standard electroweak model at accelerators. 

In fact, the standard solar model, which is only an approximate de- 
scription of the Sun, is surprisingly successful in predicting the results of 
the pioneering solar neutrino observations, in particular those obtained 
after 1986. It is our considered opinion that both the standard solar 
model and the solar neutrino experimental results consist a great tri- 
umph for theoretical and experimental physics, in spite of the so called 
"solar neutrino problem". 
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